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STABILITY ANALYSIS AND TREND STUDY OF A BALLOON TETHERED 
IN A WIND, WITH EXPERIMENTAL COMPARISONS 

By L. Tracy Redd, Samuel R. Bland, 
and Robert M. Bennett 
Langley Research Center 

SUMMARY 

A stability analysis and trend study for a balloon tethered in a steady wind are pre- 
sented. The linearized, stability-derivative type analysis includes balloon aerodynamics, 
buoyancy, mass (including apparent mass), and static forces resulting from the tether 
cable. 

The analysis has been applied to a balloon 7.64 m in length, and the results are com- 
pared with those from tow tests of this balloon. This comparison shows that the analysis 
gives reasonable predictions for the damping, frequencies, modes of motion, and stability 
boundaries exhibited by the balloon. 

A trend study for the 7.64-m balloon was made to illustrate how the stability bound- 
aries are affected by changes in individual stability parameters. The trends indicated in 
this study may also be applicable to many other tethered-balloon systems. 

INTRODUCTION 

A tethered balloon is useful for a number of purposes, such as supporting antennas 
or providing an aerial platform. Such operations are often impaired by the occurrence 
of dynamic instabilities, especially during strong wind conditions. Although information 
relating to the stability of towed and tethered bodies, including balloons, has been pub- 
lished (see refs. 1 to 4, for example), a systematic procedure for the analysis of tethered- 
balloon stability is needed. In an attempt to fill this need, the NASA Langley Research 
Center undertook a general research study to develop improved techniques for predicting 
the stability of tethered balloons. The scope of this study includes: derivation of the 
stability equations, adaptation of these equations for use in computer programs, formu- 
lation of methods for determining stability derivatives, verification of the stability analy- 
sis by the use of tow-test experiments, and an investigation of the effects of changes in 
stability parameters on the balloon stability boundaries. 



Portions of the tethered-balloon stability study are given m references 5 to 8. Ref- 
erences 5 and 6 describe the methods used to obtain the physical properties of the tow- 
test balloon. Reference 7 contains a description of the computer programs used for cal- 
culating and plotting the analytical results. A brief overview of the complete study is 
given in reference 8. 

The purpose of the present paper is to derive the equations of motion for a tethered 
balloon and to assess the accuracy of these equations by comparing the analytical results 
with results from tow-test experiments. The equations of motion include balloon aero- 
dynamics, apparent mass, structural weight, buoyancy, and static forces resulting from 
the tether cable. An additional purpose is to present a trend study which illustrates the 
effects of changes in balloon properties on the stability boundaries. 

SYMBOLS 


distance along balloon center line from nose to reference point, m (see fig. 4) 


D 


'D, 


' m 


'n 


buoyancy force, N 


drag coefficient, 


D 


PVjV 2| t 

tether -cable drag coefficient 
L 


lift coefficient, 


pVjs/2 


M 


rolling-moment coefficient, 


X,A 


pV 2 S c/2 


M- 


pitching- moment coefficient, 


Y,A 


yawing-moment coefficient, 


pVj Sc/2 
m Z,A 


pVjSe/2 


side-force coefficient, 


Y,A 


PVjS/ 2 


balloon length, m - 


2 



D aerodynamic drag force, N 

d £ tether -cable diameter, m 

^X’^Y’^Z external forces acting on balloon parallel to X-, Y-, and Z-axes, 

respectively, N (see fig. 2) 

F t,F »,F t external forces acting on balloon parallel to X* -, Y'-, and Z'-axes, 
X Y Z 

respectively, N (see fig. 1 for axis system) 

f = cos y 


h^j. component of distance from reference point to center of buoyancy, positive 

for center of buoyancy below reference point, m (see fig. 4) 


h C g component of distance from reference point to center of mass, positive for 

center of mass below reference point, m (see fig. 4) 


hi. ,h. coordinates defined by equations (A23) 

K 1 k 2 

h sr component of distance from reference point to center of mass of balloon 

structure, positive for center of mass below reference point, m (see fig. 4) 


I x ,Iy,I z rolling, pitching, and yawing moments of inertia, respectively, about balloon 
center of mass in the stability axis system (including aerodynamic apparent 
inertias), kg-m2 


Ixx.fyyj^z rolling, pitching, and yawing moments of inertia, respectively, about 

balloon center of mass in the body-reference principal axis system 
(including aerodynamic apparent inertias), kg-m2 


Ixy,Ixz>Iyz products of inertia about balloon center of mass in the XY-plane, 

XZ-plane, and YZ-plane, respectively (including aerodynamic 
apparent inertias), kg-m2 

k xx tether x-force per unit of x-displacement at bridle confluence point, N/m 

(see eq. (B21a)) 

k xz tether x-force per unit of z-displacement, N/m (see eq. (B21b)) 
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kjjg tether x-force per unit of pitch displacement, N/rad (see eq. (A27a)) 

kyy tether y-force per unit of y-displacement at bridle confluence point, N/m 

(see eq. (B24)) 

ky^ tether y-force per unit of roll displacement, N/rad (see eq. (A27h)) 

ky^, tether y-force per unit of yaw displacement, N/rad (see eq. (A27i)) 

k zx tether z-force per unit of x-displacement, N/m (see eq. (B21c)) 

k~~ tether z-force per unit of z -displacement at bridle confluence point, N/m 

(see eq. (B21d)) 

k z q tether z-force per unit of pitch displacement, N/rad (see eq. (A27b)) 

k 0x tether pitching moment per unit of x-displacement, N-m/m (see eq. (A27c)) 

tether pitching moment per unit of z-displacement, N-m/m (see eq. (A27d)) 
t/z 

k 00 total tether pitching moment per unit of pitch displacement in the body-axis 

system for pitch about center of mass, N-m/rad (see eq. (A27e)) 



portion of kgg due to rotation of balloon relative to steady tension vector 
at bridle confluence point, N-m/rad (see eq. (A27f)) 



portion of kgg due to displacement of bridle confluence point, N-m/rad 
(see eq. (A27g)) 


k^y tether rolling moment per unit of y-displacement, N-m/m (see eq. (A27j)) 

k tether rolling moment per unit of roll displacement in the body-axis system 
for roll about center of mass, N-m/rad (see eq. (A27k)) 


k , , tether rolling moment per unit of yaw displacement, N-m/rad (see eq. (A271)) 

< P V' 

k^. tether yawing moment per unit of y-displacement, N-m/m (see eq. (A27m)) 

k^ tether yawing moment per unit of roll displacement, N-m/rad (see eq. (A27n)) 


4 



k^ tether yawing moment per unit of yaw displacement in the body-axis system 

for yaw about center of mass, N-m/rad (see eq. (A27o)) 

L aerodynamic lift force, N 

l tether cable length, m 

Z br component of distance from reference point to center of buoyancy, positive 

for center of buoyancy forward of reference point, m (see fig. 4) 

l component of distance from reference point to center of mass, positive for 

eg 

center of mass forward of reference point, m (see fig. 4) 

Z sr component of distance from reference point to center of mass of balloon 

structure, positive for center of mass aft of reference point, m (see fig. 4) 

I component of distance from reference point to bridle confluence point, posi- 

tive for confluence point forward of reference point, m (see fig. 4) 

M s ,M S moments defined by equations (A31) 

1 2 

M x ,M y ,M z rolling, pitching, and yawing moments about X-, Y-, and Z-axes, 

respectively, N-m (see fig. 2) 

m mass of balloon, kg 


m a 

m ai > m a 2 > m 


m g 

m s 

m-J. 

m x ,my,m z 


apparent mass of air associated with accelerations of balloon, kg 

a apparent masses associated with accelerations in x-, y-, and 
z-directions, respectively, kg (see eqs. (A10) and (All)) 

mass of inflation gas, kg 

balloon structural mass (including bridle, test instruments, and payload), kg 
combined mass of balloon structure and inflation gas, mg + m s , kg 
masses defined by equations (A32) 
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m x a’ m y,a> m z,a aerodynamic apparent masses associated with balloon accelera- 

tion along X"-, Y"-, and Z"-axes, respectively, kg (see fig. 1 
for axis system) 

m x,o ,m y>o> m z,o total masses of balloon for accelerations in x-, y-, and 

z-directions, respectively, kg (see eqs. (A10)) 

cable drag per unit length for cable normal to the wind, N/m 

total rolling, pitching, and yawing rates about X-, Y-, and Z-axes, respec- 
tively, rad/sec (see fig. 2) 

cable drag force defined by equation (Bl) 

perturbation rolling, pitching, and yawing rates about X-, Y-, and Z-axes, 
respectively, rad/sec 

q = \jl + p 2 

S characteristic area of balloon, Vb*^, m2 

T tether -cable tension, N 

Tq.T^ tensions of tether cable at lower and upper ends, respectively, N (see fig. 1) 
t time, sec 

t^ r component of distance from reference point to bridle confluence point, posi- 

tive for confluence point below reference point, m (see fig. 4) 

U,V,W velocities of balloon center of mass along X-, Y-, and Z-axes, respectively, 
m/sec (see fig. 2) 

u,v,w perturbation velocities of balloon center of mass along X-, Y-, and Z-axes, 

respectively, m/sec 

steady wind velocity, m/sec 


P,Q,R 

Pn 

P,q,r 
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V b volume of balloon hull (i.e., gas bag), m3 

V n component of wind velocity normal to cable, sin y, m/sec 

v* y- perturbation velocity of balloon with respect to airstream, v - V^, m/sec 

W„ structural weight of balloon (including bridle, payload, and test instruments), N 

w* z- perturbation velocity of balloon with respect to airstream, w + V^d, m/sec 

w c tether-cable weight per unit length, N/m 

x,y,z perturbation displacements in the stability axis system, m (see fig. 1 for axis 

system) 

x',y',z' perturbation displacements in the earth-fixed axis system, m (see fig. 1 for 
axis system) 

x",y",z" perturbation displacements in the balloon-fixed axis system, m (see fig. 1 for 
axis system) 

I 

x c’ y c’ z c perturbation displacements of upper end of tether cable (e.g., bridle conflu- 
ence point) in the earth-fixed axis system, m 

xj distance parallel to X-axis from reference point to center of mass, positive 

for center of mass forward of reference point, m (see appendix C) 

x^,Zj coordinates of balloon center of mass with respect to tether-cable anchor 

point; xj is horizontal and positive in downstream direction, Zj is 
vertical and positive upwards, m (see fig. 1 for tether anchor point) 

z j. distance parallel to Z-axis from reference point to center of mass, positive 

for center of mass below reference point, m (see appendix C) 

a perturbation angle of attack of X-axis, ^ — + 9, rad 

OO 

balloon angle of attack at trim, rad 

y ' 

(3 angle of sideslip, ^ \p, rad 

* rv> 
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tether-cable angle with respect to the horizontal (see appendix B), rad 




angles between the horizontal and tether cable at lower and upper ends, 
respectively, rad (see fig. 1) 


displacement defined by equation (B22) 


angle between principal X-axis of balloon and stability X-axis, rad 


V 


real part of characteristic root of stability equations, damping parameter, 
1/sec 


0,<3>,'I' Euler angles of pitch, roll, and yaw, respectively, rad (see fig. 3) 
6,<p,\p perturbations of 0, $, and respectively, rad 

X characteristic root of stability equations, 17 ± ico, 1/sec 


X 

A 0’ x l 

P 


variable defined by equation (B12) 


X at lower and upper ends, respectively, of tether cable 


atmospheric density, kg/m3 


variable defined by equation (B9b) 


T 0’ T 1 


w 


t at lower and upper ends, respectively, of tether cable 

imaginary part of characteristic root of stability equations, circular fre- 
quency, rad/sec 


Subscripts: 


aerodynamic force term 


B 


buoyancy force term 


tether -cable force term 


gravity force term 
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R 


reference point (see fig. 4) 


t equilibrium trim condition 

0 lower end of tether cable 

1 upper end of tether cable 

Subscripts used with the aerodynamic coefficients denote derivatives as follows: 


a 

with respect to 

a 

a 

with respect to 

dc/2V 00 

f3 

with respect to 

(3 

P 

with respect to 

/3c/2V 00 

P 

with respect to 

pc/2V 00 

q 

with respect to 

qc/2V 00 

r 

with respect to 

rc/2V OQ 

u 

with respect to 

U /Voo 


Dots over a symbol indicate derivatives with respect to time. 

ANALYSIS 


General Comments and Assumptions 

The stability analysis presented is for a balloon tethered from an earth-fixed anchor 
point and flying in a steady wind. The balloon's motions consist of small perturbations 
about a steady-flight reference condition. 

The analysis is essentially a linearized, stability-derivative type of analysis simi- 
lar to that for a conventional airplane (ref. 9). In contrast with conventional-airplane 
analyses, however, the equations of motion for the tethered balloon include buoyancy 
forces, apparent mass terms, and static forces resulting from the tether cable. 
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The coordinate systems, forces, and moments used in deriving the equations of 
motion for the tethered balloon are shown in figures 1 to 4. The three different sets of 
right-hand Cartesian coordinates used in this analysis are as follows: 

(1) The stability axis system (X-, Y-, and Z-axes) is a balloon-fixed coordinate 
system with its origin located at the balloon center of mass. When the balloon is in the 
reference steady-state condition, the X-axis points into the direction of the wind and the 
Z-axis points downward. The XZ-plane is the plane of symmetry. 

(2) The X'-, Y' -, and Z'-axes are an earth-fixed coordinate system. When the bal- 
loon is in a reference steady- state condition, the origin of the coordinate system is located 
at the balloon center of mass. The Z'-axis points vertically downward, the Y'-axis points 
normal to the wind, and the X'-axis points forward. 

(3) X"-, Y"-, and Z"-axes are a balloon-fixed coordinate system with its origin 
located at the balloon center of mass. The X"-axis is parallel to the axis of symmetry 
of the balloon and is positive toward the nose. The x '- and Z ! '-axes are in the plane of 
symmetry, and the Y"-axis points toward the right side of the balloon. 

The mathematical model used in the analysis includes the following considerations 
and assumptions: 

(1) The equations of motion are referenced to the center of mass of the balloon. 

The balloon mass includes the mass of the balloon structure (including bridle and pay- 
load), the buoyant gas inside the balloon, and the apparent air mass; hence, the balloon 
center of mass is at a different location than the structural center of mass (see fig. 4). 

The apparent air mass is defined as the mass of air surrounding the balloon which moves 
with the balloon when it is accelerated. 

(2) There are four kinds of external forces acting on the balloon, namely, aerody- 
namic forces, tether-cable forces, buoyancy forces, and gravity forces (fig. 1). 

(3) The balloon is symmetric laterally and has yaw, roll, and sideslip angles equal 

to zero m the reference steady- state trimmed condition ^i.e., ^ ^ ^ = 0). 

(4) The steady wind velocity is parallel to the horizon. Thus for steady-state 

equilibrium conditions, the pitch angle of the X-axis is zero ^i.e., 0 t = 0^ and the stability 
and earth-fixed axes systems are coincident. 

(5) The balloon and bridle form a rigid system. 

(6) The tether cable is flexible, but inextensible, has negligible mass compared to 
that of the balloon, and contributes only static forces at the bridle confluence point. 

(7) The cable weight and drag normal to the cable are needed only for determining 
the static cable forces, equilibrium shape of the cable, and the cable derivatives. 

(8) The longitudinal and lateral equations of motion are uncoupled. 
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The preceding assumptions simplify the stability analysis, but they also restrict the 
analysis to balloons with light tether cables. One difficulty is that the tether cable effec- 
tively increases the mass of the balloon system by approximately one-third the mass of 
the cable (ref. 10). Consequently, if the effective cable mass is considered to be concen- 
trated at the bridle confluence point, then its contribution to the balloon's mass proper- 
ties can be appreciable. The mass of the tether cable can also cause coupling of the cable 
and balloon motions, but this effect is neglected in the present study. Because of these 
difficulties, it is believed that the present analysis should be restricted to cases where 
the cable mass is small compared to that of the balloon mass (including the apparent 
mass). Data which are presented later, however, indicate that the present analysis gives 
good correlation with experimental results for tether-cable masses as large as 16 per- 
cent of the balloon mass. 

Although the present analysis is somewhat restrictive, it can be used for a large 
class of tethered-balloon systems. Furthermore, the analysis yields a system of second- 
order ordinary differential equations which can be solved by methods used in conventional- 
aircraft stability analyses. Recent analyses (refs. 11 and 12) have been developed which 
include the mass and dynamics of the tether cable, but these analyses involve rather com- 
plex methods of solution. 

The equations of motion for the mathematical model used in the analysis are derived 
in appendices A and B. Note that the mass and inertia terms in these equations (see 
eq. (A10)) include the apparent mass of air associated with the balloon accelerations. In 
contrast with conventional -airplane analysis, the balloon analysis must include the appar- 
ent mass because it accounts for greater than 60 percent of the total mass for some bal- 
loon systems. 

Although the equations of motion shown in appendix A are derived for a balloon 
tethered in a steady wind to an earth-fixed anchor point, the equations are also valid for 
a balloon towed at a constant velocity over a horizontal surface. 

Equilibrium Trim Conditions 

The steady- state trimmed conditions are obtained from equations (A30a), (A30c), 
and (A30e) of appendix A by setting the perturbation quantities equal to zero. The results 
are: 


pvjs 

2 


C D - Tj cos 7l 


= 0 


( 1 ) 


PV oo 2 S 
2 


+ B - W s 


- Tj sin = 0 


( 2 ) 
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pV 2sc 

2 C m + sin - h k 2 T l cos - M s 2 = 0 


(3) 


Substituting equations (1) and (2) into equation (3) to eliminate the cable tension and 
angle y ^ gives the following trim equation: 


pVoo 2 S 


IV 2 


C L + B - w s - h kl 


^PVoo 2S o \ PVoo 2 Sc 


OO P 

4 2\ 2 D ‘ 


f — c m - m S2 = o 


( 4 ) 


where from equations (A23) and (A3 lb) 

h k x = (*tr - *cg) cos «t + ( 4 tr - h cg )sin a t 
h k 2 = (hr - h C g)cos a t - - l c g)sin a t 




Equation (4) can be solved by Newton iteration to find the equilibrium trim angle of attack 
at for various wind velocities, provided the aerodynamic coefficients C L , Cj), and 
C m are known functions of Once the trim angle of attack is determined, then equa- 
tions (1) to (3) can be solved to find and y^. The a -dependent stability coefficients 

can also be evaluated. 


Stability Equations 

The stability equations are obtained by setting the equilibrium trim portions of the 
equations of motion (eqs. (A30a) to (A30f)) equal to zero. The remaining expressions give 
six linear stability perturbation equations, three of which involve only longitudinal vari- 
ables and three of which involve only lateral variables. For convenience, the primes on 
the x', y', and z* variables are dropped. Thus, the following working forms of the 
stability equations written about the balloon center of mass are obtained. 

Longitudinal equations of motion .- The equations of longitudinal motion may be 
written as follows: 

x-force 


x + 


PVcoS 

2m y 


(2C D + C D u )i 


kxx pV S 
x + x + — “ 
m x 




z + 


^xz 


z + 


k x „ PVjsC„ 


x0 

m x 


a 


2m, 


0=0 


(5a) 
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z-force 


^£/2C l + Ct )i 
2m z { L L uJ 


k7v •• 

x + x + z + 


m. 


^fC L + C D \z + ^z Z 

2m z ( ^o! m z 


^(c L + c L 2),.o 

4m z L a L q / \ m z 2m z ) 


(5b) 


Pitching moment 


PVoqSc 


2I V \ 


(2C m + C 


\ • “0x 

)X + — — 

u / Iy 


m T 


x - 


pSc 2 .. pVooSc 


4L, 


C m .z - 


a 21y 


K fl7 

C m Z + -T — Z + 0 

m a I v 


PV^SC 2 


'M 


si pV^Sc 


( Cm (i + Cm q)®' + 1 ' Tf 2Iv ^ ' I, 


c mQ! + tT = 0 
■y aL y 


(5c) 


Lateral equations of motion .- The equations of lateral motion may be written as 
follows: 


y-force 


y - 


pVooS kyy pv^sc 


2m, 


c Y«y + m y - 
y P m y 

P Y //-> r> \ ", 

T^( CY fi ' N* + 


C y + 


4m y ~ Y p 


fa <P pV 0O 2 SC L 


m, 


2m, 




' P vjs(c Yp + c D ) kyi ; 


2m, 


m. 


\jy = 0 


(6a) 


Rolling moment 

PSc 2 PVqqSc 

4I X ^ 2I X 


<£y 


V + 1 T y + < ^ 


PV ooSc 2 


4L 


C, 4> 


k 4>4> * Vi sln n + “*! , I xz v 

f <P - T~ W 

lx lx 


PV^Sc 2 


+ T(V C ^' + 


2I V 


\f/ = 0 


(6b) 


13 



Yawing moment 


pSc 2 » PVcoSc ■ Vy Ixz v pVcqSc 2 c i 
n /3 y 2I Z % Y + Iz y Iz 0 4I Z 


k i//0 + M s, ' h k 1 T l sin .. pV^Sc 2 , 

— V 1 * ♦ * + T(% • c "r * 


+ Vd) 


21 , 


\p = 0 


(6c) 


Although the stability equations and equilibrium trim equations are written about 
the balloon center of mass, the aerodynamic forces and moments may be referenced to 
an arbitrary point (see the reference point in fig. 4) which is generally not located at the 
center of mass. In such cases it is necessary to transfer the aerodynamic terms to the 
balloon center of mass through the use of the transfer equations given in appendix C. 

The aerodynamic forces are treated in this manner in the present study because it 
is more convenient to use aerodynamic inputs to the computer programs (ref. 7) which 
are referenced to a fixed point in the balloon rather than to a variable center of mass. 
The location of the balloon center of mass varies with changes in altitude because the 
mass of the buoyancy gas and apparent mass are both functions of atmospheric density. 
The location of the center of mass also changes as the payload position and weight are 
varied. 


Balloon Stability Characteristics 

The longitudinal and lateral stability characteristics of the balloon are determined 
independently from equations (5) and (6), respectively, in a manner similar to that used 
in conventional-airplane stability analysis (ref. 9). Since the stability equations are ordi- 
nary linear differential equations with constant coefficients, they have transient solutions 
which are always exponential in form. For example, a typical variable such as 9 is of 
the form 9 = 9e^, where 9 is a complex constant. If these exponential forms are sub- 
stituted into the longitudinal and lateral stability equations and the determinants of the 
coefficients are set equal to zero, solutions for the characteristic roots X can be 
obtained. 

Solving equations (5) and (6) in this manner gives six characteristic roots for each 
of the longitudinal and lateral cases. These roots appear as complex conjugate pairs 
(i.e., X = 77 ± icu) for oscillatory modes of motion or as real numbers (i.e., X = 77) for 
the aperiodic modes. Thus, in each of the longitudinal and lateral cases, it is possible 
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for the balloon to exhibit from three to six modes of motion, depending on whether the 
roots are real or complex conjugate pairs. For any mode, the balloon system is stable 
when Tj < 0, neutrally stable when 77 = 0, and unstable when r\ > 0. 

Computer programs based on equations (5) and (6) were developed for calculating 
the stability characteristics and plotting the results. These programs are described in 
detail in reference 7. Some typical results from these programs are given in the section 
entitled "Results and Discussion." 


TOW TESTS 

In order to provide experimental data for comparison with analytical results, a bal- 
loon was towed from a truck. Some instabilities were measured and some frequencies 
and dampings were obtained for the stable cases. 

Description of the Balloon and Tow- Test Apparatus 

Balloon .- The balloon used for the tow tests is shown attached to the tow truck in 
figure 5. The balloon is 7.64 m long and has a basic shape similar to the Navy C-class 
airship configuration (ref. 13), except that the aft section is conical and the nose is nearly 
spherical. The essential components of the balloon consist of an inflatable hull and a 
rigid (i.e., not inflated) tail-fin assembly. A load band is attached to the lower portion 
of the hull which permits variations in the bridle attachment points and, hence, makes it 
possible to change the balloon equilibrium trim angle of attack. A more detailed descrip- 
tion of the balloon is given in reference 6. 

The aerostatic, mass, inertia and aerodynamic properties of the balloon for the 
reference configuration are given in table I. The reference configuration used in the 
present study was arbitrarily chosen such that the balloon would fly stably at small trim 
angles of attack for a large range of tow velocities. The balloon properties shown in 
table I were determined by procedures described in reference 6. 

Test equipment and procedures .- The panel truck used in the tests (fig. 5) provided 
space for data recording equipment and acted as a towing vehicle. The test instrumenta- 
tion on the truck included a load cell to measure cable tension, potentiometers to measure 
cable inclination angles, and an anemometer (i.e., wind sensor) to measure speed and 
sideslip angle. Instruments were also suspended rigidly beneath the balloon to measure 
pitch and roll motions. Outputs from the instruments were recorded on strip charts 
during the tests. Motions of the balloon were also recorded by a camera mounted on 
the truck. 

There were several features of the test equipment and procedures which limited the 
range of the tow tests. The primary limitation was that 61m was as long a tether cable 
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as was practical or convenient to use. The maximum speed of the tow truck with the bal- 
loon attached was 29 m/sec. 


Test Description and Data Reduction' 

A list of the tow tests and test variables is given in table II. The test variables are 
tow velocity V^, bridle-confluence-point location ^ r , and tether-cable length l, weight 
w c , and diameter d c . The balloon trim angle of attack depends upon the tow veloc- 
ity and bridle- confluence-point location. The range over which these properties are 
varied is indicated in table n. 

Each test result shown in table II represents the average of several runs at essen- 
tially the same relative airspeed. These runs were made during periods of calm air when 
the ambient wind velocity was less than 2.5 m/sec. 

During tests 10 to 14, the balloon was initially displaced with an auxiliary line from 
its equilibrium trim angle and released when the tow vehicle was at the proper speed. 

This was done in order to produce disturbances from which the frequency and damping 
of the balloon motion could be determined. The oscillations observed in tests 15 to 18 
were self-excited unstable motions requiring no initial displacements. 

The pitch angle, roll angle, velocity V w , and tether-cable data were recorded on 
the strip charts. The x-, y-, and (//-displacements of the balloon were obtained from the 
movie film. No attempt was made to measure the z-displacements. 

RESULTS AND DISCUSSION 

Since the longitudinal and lateral equations of motion are uncoupled, the longitudinal 
and lateral results are presented separately. In each case, the results for the 7.64-m 
balloon in the reference configuration (i.e., the configuration for tests 12 and 15 of table II) 
are presented first. These results are followed by data for the balloon in the other test 
configurations (i.e., the configurations for tests 10, 11, 13, 14, 16, 17, and 18). A param- 
eter trend study is also presented for both the longitudinal and lateral cases. 

Where possible, the results from the analysis are compared with experimental data 
from the tow tests. The quantities compared are the oscillatory frequencies, damping 
rates, modes of motion, and stability boundaries for the 7.64-m balloon. 

Longitudinal Stability Characteristics 

Reference configuration .- The calculated values of the longitudinal frequencies co 
and damping rates 77 for the balloon in the reference configuration are plotted as a func- 
tion of wind velocity in figure 6(a). These same values of u> and 77 are given m root- 
locus form with velocity as the varying parameter in figure 6(b). These figures indicate 
that the reference configuration has three oscillatory modes of motion at velocities less 
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than 13.5 m/sec. For velocities greater than this value, mode 1 splits into two real non- 
oscillatory modes. The figures also show that the calculated value for mode 2 becomes 
unstable (i.e., rj = 0) at 22.4 m/sec. 

The measured results from the tow tests of the balloon in the reference configura- 
tion (i.e., tests 12 and 15 of table II) are shown in figure 6(a). In comparing these results 
with the calculated data, it appears that the measured frequencies and damping rates cor- 
respond to those calculated for mode 2, except at a velocity of 12.5 m/sec. At this veloc- 
ity, two different frequencies were observed - one corresponding to mode 2 and the other 
to mode 3 - but the data were insufficient to determine the damping rates for either mode. 

The measured data show the balloon to be unstable at velocities greater than 
24 m/sec, which is in good agreement with the calculated results. The instabilities 
observed at these velocities were quite violent. In several cases the motions become 
so excessive that either the tether cable or load band broke and the balloon was severely 
damaged. 

As a means of illustrating the longitudinal modes of motion, the tethered balloon and 
its center of mass are shown for a sequence of time intervals in figures 7(a), (b), and (c). 
The plots in these figures are drawn from the calculated data by means of a computer 
program given in reference 7. Except for the middle and lower portions of figure 7(a), 
the time sequence shown for each velocity and value of X represents one cyle of free 
oscillation for a single mode. The time sequences for the plots in the middle and lower 
portions of figure 7(a) are just long enough to allow the nonoscillatory balloon motions to 
reach steady- state equilibrium conditions. 

Each mode of motion shown in figure 7 was obtained by substituting the value of the 
characteristic root X of the mode into equation (5) and solving for the ratios of the dis- 
placements to the pitch angle. The initial pitch-angle amplitude was an input to the com- 
puter program (ref. 7) and its value could be arbitrarily selected. Thus the initial ampli- 
tude of the balloon's displacements could be adjusted as desired. 

The modes of motion shown in figure 7 do not represent transient response in the 
usual sense. During transient motions, all of the modes of a dynamic system generally 
participate in the motion simultaneously, the relative degree of participation depending 
on the initial conditions or excitation. The purpose of figure 7 is simply to illustrate the 
character of each mode individually. When one of the modes becomes lightly damped or 
unstable, however, it may dominate the transient motions. Such was the case in the runs 
of test 15 for velocities greater than about 20 m/sec. At these velocities the observed 
transient motions of the balloon were very similar to those shown for mode 2 in the lower 
portion of figure 7(b). 

The effect of damping on a given mode, at a particular velocity, can be observed in 
figure 7 by comparing the balloon displacements at the beginning and end of each time 
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sequence. For example, in figure 7(b) it can be seen that the balloon motion is highly 
damped (i.e., r/ < 0) at a velocity of 10 m/sec since the displacement is almost zero at 
the end of the time sequence. At = 30 m/sec, it is observed that the balloon is unsta- 
ble (i.e., rj > 0) since the displacement at the end of the sequence is larger than that at 
the beginning of the cycle. 

It was mentioned previously that mode 1 splits into two nonoscillatory modes (i.e., 
modes la and lb) at velocities above 13 m/sec. Figure 7(a) shows that mode la is pre- 
dominantly translational motion while mode lb is predominantly rotational motion. Fig- 
ure 7(c) indicates graphically the result presented in figure 6(a), that is, the balloon 
motions for mode 3 become more damped as velocity increases. 

A computer generated motion-picture supplement (L-1118) has been made which 
shows the balloon modes of motion in more detail for both the longitudinal and lateral 
cases. This film is available on loan and a request card is included at the back of this 
report. 

Additional configurations .- Results for the configurations corresponding to tests 10, 
11, 13, 14, 16, 17, and 18 of table II are presented in figures 8 to 13. These figures show 
that the calculated values of oj and rj are in fair agreement with the measured values 
at the test velocities indicated. Note that in some of the figures there are fewer measured 
values of r] presented than of w. For those cases where rj is omitted, the measured 
data were inadequate to determine the damping of the balloon motions. 

The primary purpose of tests 10 to 14 was to check the accuracy of the frequencies 
and damping rates calculated from the analysis. An additional objective was to determine 
experimentally the effects of the tether-cable length, diameter, and weight on the modes 
of motion. Neither of these test objectives required the use of high tow velocities; thus, 
tests 10 to 14 were conducted at velocities well below the predicted velocities of instability 
(see figs. 6 and 8 to 11). 

By comparing the data from tests 10, 11, and 12 (figs. 8, 9, and 6 and table n) it can 
be seen that changes in the tether-cable length have large effects on all three modes of 
motion. These data show that the velocity at which mode 1 becomes nonoscillatory 
decreases as the tether length increases. The velocity at which mode 2 becomes unsta- 
ble also decreases as the tether length increases. (The trend study for the present bal- 
loon configuration, which is presented later, indicates that the velocity of instability 
increases as the tether length increases for lengths longer than 90 m.) At a given veloc- 
ity, the frequencies of all the modes decrease as the tether length increases, but mode 3 
shows the greatest change in frequency. 

The effects of changes in the tether-cable weight and diameter on the longitudinal 
modes of motion can be determined by comparing the data for tests 11, 13, and 14 (figs. 9 
to 11 and table II). These data show that the velocity at which mode 2 becomes unstable 
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decreases as the cable weight and diameter increase. Based on the trend study, which is 
presented later, the modes of motion are affected more by changes in the tether-cable 
diameter than by changes in cable weight. It should be noted that the tether cable for 
tests 14 and 18 (fig. 11) was the heaviest cable used in the tow tests. The mass of this 
cable was about 16 percent of the balloon mass (including apparent mass). 

The data presented in figures 11 to 13 (and also fig. 6) indicate that the stability 
analysis does predict the damping and frequencies near the point of instability (i.e., near 
7j = 0) quite well. Note that in figures 11 and 13 the slopes of the damping curve are 
almost flat near r; = 0. For these two cases a slight error in r? makes a large differ- 
ence in the -velocity at which the balloon goes unstable. This result is particularly true 
for figure 13. For this case the calculated and measured damping values seem to corre- 
spond quite well. The calculations, however, indicate that the motions are very lightly 
damped and may never go unstable, whereas the measured results show that the motions 
become slightly unstable for velocities greater than about 23 m/sec. It is possible for 
this case that the measured damping may be slightly in error since it was difficult to 
determine from the recorded motions whether the balloon was very lightly damped, neu- 
trally damped, or slightly unstable. There is some evidence to indicate that the motions 
may have been very lightly damped since the cable loads during the test did not become 
excessive and the motions were quite mild compared to those in test 15. 

The calculations shown in figure 12 are for the case where the balloon is tethered 
at a high angle of attack (see table II). For this case the balloon is unstable at low wind 
velocities. This unstable condition is consistent with the results of reference 11 and is 
primarily due to interactions between the balloon's lift and drag forces. For this condi- 
tion the lift acts as a forcing function and the drag acts as a damping function, and thus a 
high lift-drag ratio aggravates the balloon's motions (ref. 11). The two experimental 
points shown in figure 12 indicate that the balloon was not quite unstable as analytically 
predicted, but was very lightly damped. The motions were quite large, however, and may 
have been either lightly damped or limit- cycle oscillations which are considered unstable 
motions. The analysis predicted that the balloon would be more stable at higher veloci- 
ties because, for this case, the angle of attack decreases with increasing velocity. Never- 
theless, because the balloon was experiencing quite large motions at such low velocities 
and visual observations gave evidence that this configuration may be unstable, no test 
velocities higher than those shown were attempted. 

Longitudinal trend study .- Analytical trend studies for the 7.64-m balloon in the ref- 
erence configuration are presented in figures 14 to 27. These figures show how changes 
in individual stability parameters affect the frequency and velocity of instability. The 
properties of the test balloon in the reference configuration are indicated in each of the 
figures. 
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The parameters considered in the longitudinal- stability trend study, and a list of 
the corresponding figures in which these parameters appear, are given in table ID. A 
preliminary investigation indicated that only those parameters shown have a major effect 
on the longitudinal stability of the balloon in the reference configuration. Notice that 
several of the aerodynamic parameter symbols shown in the table have the prefix A. 

This prefix indicates an incremental difference between two values of the same parame- 
ter, where one of the values is for the balloon in the reference configuration ^e.g., 

ACd - Cj) - C-D re f erence )* 

The effects of the balloon mass and inertia properties on the stability boundaries are 
shown in figures 14 to 17. These figures indicate that a reduction in the balloon mass 
m«j> and pitch inertia should raise the instability speed. They also show an increase in 
stability speed (fig. 16) with the center of mass located lower and farther forward on the 
hull. Further, figure 17 indicates that reducing m x a and increasing m z>a would 
increase stability. 

The effects of buoyancy are shown in figure 18. This figure indicates that for better 
stability the buoyancy force should be larger and located farther forward in the hull, but 
not so far forward as to get into low speed instabilities. 

The effect of the tether-cable length on the stability boundary has been discussed 
previously in relation to figures 6, 8, and 9, and additional information is presented in 
figure 19(a). A plot showing the equilibrium position of the balloon as a function of tether 
length and wind velocity is given in figure 19(b). Note that in figure 19(a) the measured 
points of instability for test 15 are compared with the calculated data. Because of the 
previously mentioned limitations in the test equipment, points of instability were mea- 
sured only for the 61-m tether length. A number of tow tests were conducted using 
tether cables of shorter lengths; but because the maximum speed of the tow truck was 
29 m/sec, no additional longitudinal instabilities were observed for the balloon using the 
same bridle arrangement and cable diameter as the reference configuration. 

Figure 20 shows that the location of the tether- cable attachment point (i.e., bridle 
confluence point) has a large effect on the stability boundary and balloon trim angle of 
attack. The analytical results in this figure show that the instability speed increases as 
the attachment point is moved forward. In contrast with the analytical results, the mea- 
sured data for tests 15 and 17 indicate that the instability speed may decrease slightly 
with increasing ^ r - As discussed earlier, however, the measured instabilities shown 
for test 17 may be somewhat questionable since it is possible that the balloon motions 
may have been very lightly damped for the velocities shown. Note that in figure 20(b) 
the trim angle decreases with increasing velocity for small values of l^ v , but the oppo- 
site is true for large values of ^ r . Thus, ideally, there is a value for l^ T for which 
the trim angle is independent of wind speed. 
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The effects of tether -cable diameter and weight per unit length are shown in fig- 
ures 21 and 22. The cable diameter affects the cable drag which, in turn, has a large 
effect on the stability boundaries, as shown in figure 21. A comparison of figures 21 and 
22 indicates that the cable diameter has a much greater effect on the stability boundaries 
than the cable weight. These figures also show that the tether cable should be small in 
diameter and light in weight. 

Plots showing the variation of the balloon's position with changes in wind velocity 
for several different cable diameters are shown in figures 21(b) and (c). These plots 
indicate that for velocities greater than 35 m/sec, the position of the balloon does not 
change appreciably as the wind velocity increases for the cable diameters shown. On 
the other hand, the balloon position does change considerably for high wind velocities as 
the cable diameter changes. 

The results shown in figure 23 indicate that the longitudinal stability boundary is 
sensitive to variations in the tether-cable derivatives. Thus, these derivatives must be 
calculated accurately. Note from appendix B that these derivatives are not independent 
parameters, but are functions of the wind speed, balloon aerodynamics, buoyancy forces, 
and tether-cable properties. Thus, since these derivatives are not independent parame- 
ters, the absolute values of the derivatives are not shown. A more meaningful interpre- 
tation of the sensitivity of these derivatives is obtained by plotting them in terms of a 
fractional variation from their normal values for the balloon in the reference condition. 
The normal reference values of the tether-cable derivatives change, of course, as the 
wind speed changes. 

The effects of the aerodynamic terms on the longitudinal stability boundaries are 
shown in figures 24 to 27. Since the reference values of the first four terms are func- 
tions of the wind speed (see table I), the absolute values of these terms are not shown, 
but incremental changes with respect to their reference values are plotted. The fifth 
term C mo! is also presented in terms of incremental changes. Figure 24 indicates 
that positive increments of static aerodynamic drag and pitching moments about the cen- 
ter of mass are stabilizing. A reduction in lift should also be stabilizing. Although a 
large increase m lift appears to be stabilizing, a large lift coefficient may cause struc- 
tural problems at high speeds. Reduction in the values of CL a and C mQ! (fig- 25) 

leads to an increase in the speed of instability. The derivative C m . (fig. 26) should 

Q!,R 

be less negative and C m (fig. 27) should be more negative for increased stability. 

Q ,R 

Note that these latter two pitching-moment derivatives are given about the reference 
point (fig. 4) instead of about the center of mass. 

The preceding aerodynamic data indicate that for better stability the balloon's hori- 
zontal tail fin should be located farther aft of the center of mass. This configuration will 
give a more stable combination of values for Cl , ^m a > C m • > anc * C m ^. The hori- 
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zontal fin should also have a small negative angle of attack with respect to the body center 
line in order to decrease lift and give a more positive pitching moment. 

Lateral Stability Characteristics 

Reference configuration .- The calculated values of the lateral frequencies co and 
damping rates 77 for the balloon in the reference configuration are plotted as a function 
of wind velocity in figure 28(a), and in root locus form in figure 28(b). The data in these 
figures indicate that none of the modes for this configuration are unstable, but mode 2 is 
very lightly damped. This result agrees with the tow-test data shown in figure 28(a). 

Note that mode 1 splits into two nonoscillatory modes of motion at a velocity of 2 m/sec 
and that mode 3 has a rather high oscillatory frequency. 

Time sequential drawings of the balloon and its center of mass for calculated lat- 
eral modes of motion are shown in figure 29. This figure is plotted in a similar fashion 
to that of figure 7 for the longitudinal analysis, except that the bottom view of the balloon 
is drawn and the magnitudes of the displacements are proportional to the yaw angle \p. 
The plots are drawn by the computer program documented in reference 7 and include the 
effects of balloon rolling motions and trim pitch angle. Notice in figure 29(a) that mode 1 
has split into two nonoscillatory modes of motion for the velocities shown. 

Additional configurations .- Results for the configurations corresponding to tests 10, 
11, 13, 14, 16, and 17 of table II are presented in figures 30 to 35. These figures show 
that the measured data correspond to the calculated results for mode 2. No experimen- 
tal points were measured corresponding to modes 1 or 3 because mode 2 dominated the 
transient motions of the balloon. 

The calculations for mode 2 of test 10 (fig. 30) show a small region where the bal- 
loon should have been slightly unstable. Although there were no experimental lateral 
instabilities observed, the measured data in figure 30 show that the motions were very 
lightly damped. Furthermore, these motions were observed to become quite large at 
times and could have been either lightly damped or limit- cycle oscillations. 

The effects of changes in the tether-cable length on the lateral modes of motion can 
be seen by comparing the data from tests 10, 11, and 12 (figs. 30, 31, and 28 and table II). 
Except for mode 2, the data show that these effects on the lateral modes are similar to, 
but less pronounced than, those discussed earlier for the longitudinal case. For mode 2, 
an increase in tether length does not decrease the velocity at which the balloon becomes 
unstable as in the longitudinal case but causes the damping of the mode to increase at 
lower velocities and decrease at higher velocities. 

Changes in the tether- cable weight and diameter have very little effect on the lat- 
eral modes of motion, as shown by the comparison of the data from tests 11, 13, and 14 
(figs. 31, 32, and 33 and table II). Note that there are no values of co and 77 shown 
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in figure 33 for test 18 because the balloon went unstable longitudinally for this test and 
no specific lateral modes of motion were experimentally detected. 

The measured damping shown in figure 34 for test 16 is less than analytically pre- 
dicted. This result is possibly due to coupling between the longitudinal and lateral modes 
of motion for this configuration. 

Lateral trend studies .- Analytical trend studies for the 7.64-m balloon are presented 
in figures 36 to 44. These figures are presented in a similar format to those for the lon- 
gitudinal trend studies. A listing of the parameters considered in the lateral trend study 
is shown in table IV. Except for the term kyy, which is not shown, only those parame- 
ters presented in the table have an appreciable effect on the lateral stability boundaries 
of the balloon m the reference configuration. Nonzero values of kyy have little effect 
on the lateral stability, but at kyy = 0 the balloon is unstable at all velocities. This 
result indicates that there must be a lateral restraining force provided by the tether 
cable even though this force can be very small. 

The effects of the balloon mass and inertia properties on the lateral stability bound- 
aries are shown in figures 36 to 38. These figures indicate that the stability can be 
improved if the value of the yaw inertia and apparent mass my )a are reduced. The 
center of mass should also be moved farther forward for increased stability. 

The effects of buoyancy are shown in figure 39. Figure 39(a) shows that the balloon 
is only unstable in a very small region at low velocities. Figure 39(b) shows that the 
present balloon is stable for all values of £ br greater than about 0.9. 

r 

The analytical results presented in figure 40 indicate that the balloon has a small 
area of instability for short tether lengths. As mentioned previously, however, no insta- 
bilities were found, but large motions were encountered which could indicate the possi- 
bility of limit-cycle oscillations or low damping for the short tether-cable conditions. 

The effects of the static aerodynamic terms on the lateral -stability boundaries are 
shown in figure 41. These results are similar to the longitudinal case in that increasing 
the drag and decreasing the lift cause an increase in the balloon's speed of instability. 

The effects of ACy^, AC n ^, and AC nr on the lateral stability boundaries are 
shown m figures 42 to 44. In general, the terms Cy^, C n ^, and Cn r must be inher- 
ently stable ^i.e., Cy^ and C nr must be negative and C n ^ must be positive j if the 
tethered balloon is to have good lateral stability. Figure 43 shows, however, that C n ^ 
can be slightly negative and the balloon will still be marginally stable because of the 
restoring moment from the tether cable. The results of figures 42 to 44 indicate that 
for increased stability the balloon should have a larger vertical tail fin located farther 
aft of the center of mass. 
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The results presented in figure 43 show that the roll degree of freedom 0 has little 
effect on the stability boundary generated by changes in C n ^- Thus, it appears that the 
roll degree of freedom has only a moderate effect on the stability characteristics of the 
present balloon configuration. 


CONCLUDING REMARKS 

A stability analysis and trend study for a balloon tethered in a steady wind have 
been presented. Equations of motion are derived which include balloon aerodynamics, 
aerodynamic apparent mass, structural mass, buoyancy, and static forces resulting 
from the tether cable. The analysis was applied to a balloon 7.64 m in length and the 
results were compared with stability data from tow-test experiments. This comparison 
showed that the analysis gave reasonable predictions for the damping, frequencies, modes 
of motion, and stability boundaries exhibited by the balloon. Based on these results, it is 
concluded that the equations and techniques presented in the investigation offer a system- 
atic and useful method for the analysis of the stability of many tethered-balloon systems, 
particularly those with relatively light tether cables. 

The trend studies presented were made by varying individual parameters for the 
7.64-m balloon. Significant effects of the magnitude and location of the buoyancy force 
and the aerodynamic apparent masses were noted. For the longitudinal case, stabilizing 
effects were also noted as a result of reducing the tether-cable diameter (or drag coeffi- 
cient), using a more negative C mQ , (derivative of pitching moment with respect to angle 
of attack), increasing the balloon drag coefficient, and decreasing the aerostatic mass. 

As the tether-cable length was varied, the minimum longitudinal instability speed occurred 
at an intermediate length, with higher speeds of instability for longer or shorter cable 
lengths. Fore and aft location of the tether attachment point also had large effects on 
both trim angle of attack and longitudinal instability. For the lateral case, stable direc- 
tional stability was necessary for good stability characteristics. For short tether lengths, 
a small region of lateral instability was calculated for low velocities. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., May 23, 1973. 


24 



APPENDIX A 


DERIVATION OF THE EQUATIONS OF MOTION 
Rigid- Body Equations of Motion 

The kinematic and dynamic equations for rigid-body motion, given in chapter 4 of 
reference 9, are used in deriving the equations of motion for the tethered-balloon system. 
The equations of reference 9 are derived about the center of mass in the stability axis 


system and are as follows: 



Summation of forces 



EF X = m(u + QW - RV) 


(Ala) 

EF Y = m(v + RU - PW) 


(Alb) 

ZF Z = m(\V + PV - QU) 


(Ale) 


Summation of moments 

£M X = I X P - I XZ R + QR(I Z | - I y ) - I XZ PQ 

SMy = Iy$ + RP(l x - I z ) +,, I xz (p 2 - R 2 ) 

SM Z = -l xz £ + I Z R + -PQ^fy ~ lx) + fxzQ^ 

where Ij^y = Iy Z = 0 because of symmetry. 

Angular velocities 

P = <j> - \j> sin O 

Q = © cos $ + ^ cos 0 sin 4> 

R = cos © cos 4> - 0 sin $ 

Linear velocities 

= U cos 0 cos 'I' + V(sin 4> sin 0 cos ^ - cos $ sin 
dt 

+ W(cos 4> sin 0 cos ^ + sin 4> sin 4?) (A4a) 


(A 2a) 
(A2b) 
(A2c) 

(A3a) 
(A 3 b) 
(A 3 c) 
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APPENDIX A - Continued 

= U cos 0 sin + V(sin sin 0 sin <!? + cos $ cos 
dt 

+ W(cos <f> sin 0 sin - sin $ cos 4-0 (A4b) 

= -U sin 0 + V sin 4? cos 0 + W cos $ cos © (A4c) 

dt 

The axis systems, motions, forces, and moments included in the above equations are 
shown in figures 1, 2, and 3. 

The preceding equations of motion can be adapted to describe the motions of a 
tethered balloon by using the assumptions given in the section entitled "Analysis" of this 
paper. These assumptions require that the balloon center of mass has no displacements 
or velocities for equilibrium trim conditions. Thus, the only movements which the bal- 
loon exhibits are small perturbations about its initial trimmed position. Hence, for small 
displacements, the following linearized relationships are obtained (ref. 9): 


U = u; V = v; 

W = w 

(A 5a) 

0 = 0; 

* = ft 

<f> = <p 

(A 5b) 

P = p; 

Q = q; 

R = r 

(A 5c) 


where the lower-case letters denote small perturbation velocities and displacements of 
the center of mass. 

In contrast with conventional-airplane analyses (ref. 9), the expression for the 
x-component of velocity U in equation (A5a) does not include the equilibrium velocity 
term V . The reason for the difference is that in the conventional analysis the airplane 
moves through still air at a speed equal to the equilibrium velocity plus the perturbation 
velocity, whereas in the present analysis the balloon is tethered from an earth-fixed 
anchor point and the only movements of the balloon are small perturbation motions. The 
airstream does flow over the balloon at a constant equilibrium velocity in the present 
analysis, but the airstream velocity only affects the magnitude of the aerodynamic forces 
and is not an explicit part of the motion of the center of mass. 

Using equations (A5) in equations (A3) and (A4) and retaining only first-order per- 
turbation terms gives 

P = <f>; q = 6; r = p (A6) 

x' = u; y' = v; z' = w (A7) 
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APPENDIX A - Continued 

Substituting equations (A5), (A6), and (A7) into equations (Al) and (A2), and retaining only 
first-order terms, results in the following linear equations of motion: 


EF X = mx' 


(A8a) 

M 

1 

*-< 

II 

3 

«<: 


(A 8b) 

2F Z = mz' 


(A8c) 

2M X = l x 4> - I xz ip 


(A 9a) 

2 My = Iy0 

y 

(A 9b) 

SM Z = "Ixz^ + 


(A 9c) 


Apparent masses and moments of inertia .- The balloon system mass term m used 
in equations (A8) includes the structural mass of the balloon ms, the mass of the gas 
inside the balloon mg, and the apparent mass m a of the air associated with the accel- 
erations of the balloon. Since the apparent mass of streamlined bodies and prolate spher- 
oids accelerating in the axial direction is different from that for accelerations in a lat- 
eral direction (ref. 14), the following relationships for the present balloon system are 
introduced: 


m x , 0 = m s + mg + m a ^ (AlOa) 

m y>0 = m s + m g + m a2 (AlOb) 


m 


z,o 


= m s + 


m„ + rm 


(AlOc) 


where m x Q , my 0 , and m z Q are the total masses and m a ^, m^, and m a ^ are 
the apparent masses associated with accelerations in the x-, y-, and z-directions, 
respectively. 

The apparent mass of the balloon moving in the x- and z-directions also depends 
upon the equilibrium trim angle of attack such that 


m a 1 = m x,a 
m a 2 = m y,a 
m a 3 = m x,a 


cos2q!j. + m z>a sin2o!f. 


sin2a!j. + m z a cos2q:^ 


(Alla) 
(A lib) 
(A 11c) 
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where m x a , my a , and m z a are the apparent masses of the balloon accelerating 
along the X”-, Y"-, and Z"-axes. 

Because the apparent mass depends upon the direction of the acceleration, the 
expressions given in equations (A8) are written 


2F X = m x>0 x’ 

(A12a) 

SFy = m y ( oy' 

(A 12b) 

EF Z = m Z)0 2' 

(A 12c) 


where the expressions for m x 0 , niy )0 , and m z o are given in equations (A 10). 

The mass moments of inertia given in equations (A9) also include the inertias of 
the balloon structure, the gas inside the balloon, and the apparent mass of the surround- 
ing air. Furthermore, these moments of inertia depend upon the orientation of the bal- 
loon as described by the following equations: 


I x = Ixx cos2e + i zz sin2 e 
Ixz = j^xx * Iz Z )sin 2e 
J y = % 

Iz = Ixx sin2e + Izz cos2 e 


(A 13a) 
(A13b) 
(A13c) 
(A 13d) 


where I^, Iyy, and I zz are the mass moments of inertia about the principal axes and 
e is the angle between the principal X-axis and the stability X-axis. In the present anal- 
ysis, the X"-, Y"-, and Z"-axes are considered to be the principal axes; hence, e = a^. 

Forces and moments acting on the balloon .- There are four different sources of 
external forces and moments acting on the balloon as shown in figures 1, 2, and 4. These 
are: (1) aerodynamic forces, (2) tether-cable forces, (3) buoyancy forces, and (4) gravity 
forces. Therefore, the equations of motion for the balloon as given in equations (A12) and 
(A9) can be written as 


X,A + F X,C + F X,B + F X,G = m x,o*’ 

(A 14a) 

Y,A + f y,c + f y,b + F Y,G = m y,oy' 

(A 14b) 

Z,A + F Z,C + F Z,B + F Z,G = m z,o'z' 

(A14c) 
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m X,A + M X,C + m X,B + m X,G = l x4> - hiz^ (A 15a) 

My ^ + My q + My g + My q = Iy0 (A15b) 

^Z,A + ^Z,C + ^Z,B + ^Z,G = '^xz 4> + (A 15c) 

where the A, C, B, and G subscripts denote aerodynamic, tether-cable, buoyancy, 
and gravity terms, respectively. The expressions for each of the external forces and 
moments are derived in the following four sections. 

Aerodynamic Forces and Moments 

The aerodynamic forces and moments about the balloon center of mass are derived 
in essentially the same manner as used in reference 9. In deriving these aerodynamic 
terms, the following assumptions are made: 

(1) The aerodynamic forces and moments are functions of the instantaneous val- 
ues of the perturbation velocities and accelerations of the balloon with respect to the 
airstream. 

(2) The perturbation velocities of the balloon with respect to the airstream are 
u, v*, and w*. 

(3) The longitudinal and lateral forces and moments are independent. 

(4) All aerodynamic derivatives with respect to acceleration quantities are 
neglected, except F z ^ + , My^*, F y ^, 1^*, and M z ^. 

(5) The F y derivative is very small and may be neglected (ref. 9). 

(6) The equilibrium trim values of the longitudinal aerodynamic forces are 

F X, A t = -D t and F Z,A t = - L t- 

(7) The equilibrium trim values of the lateral aerodynamic forces and moments 

F Y, A t ’ M X,A t ’ and M Z,A t are zero. 

Using the preceding assumptions and equation (A6), the aerodynamic forces and 
moments are: 


X,A = - D t + F X U U + F X W * W * 

(A 16a; 

Y,A = F Y v * v * + + F Y/ + F Y 

(A 16b) 

'■7 A = -L t + F 7 u + F 7 w* + F 7 w* + F 7 0 

*‘> A 1 ^u z w* ^q 

(A 16c) 
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where 


m x,a = M x v * v * + M x^/* + M Xp^ + M x/ 

My ^ = My ^ + M Y u + ^Y w * W * + M Y^ * + My^0 


M 


Z,A = M Z y * v * + M Z- j* + M Z p 4> + M Z r ^ 


3F 


X,A 


9M 


X u 3u 


, M v . = 


Y,A 


W * 3w 1 


, My 


9M y,a m _ 9M z,a 


3q ’ 


z - 


3p ’ 


(A16d) 

(A16e) 

(A16f) 


The first-order approximations of the perturbation angle of attack a and sideslip 
/3 are (ref. 9): 


a = 


w 

V, 


v* 

^ _ v - 

V OO 


Applying this approximation to equations (A16) and expanding a in terms of z'/V^ 
and 9, and /3 in terms of and ij/, gives the following results: 

F-y- a = -D t + F y i' + Fy. + Fy- 9 (A 17a) 

a, a L a u A a Voo 


F Y,A Fy ^ + F Y j3 V w 


+ Fy p 0 + ^Fy r - Fy^lj/ - Fy^ 


F Z.A = - L t + F Z U *' + F Z i ^ + F z a ^ + (% + Fz a) (i + F V 

Mv a = My . ~ — + My ~ — + My 4) + /M y - My - My ip 
X,A P V co X /3 V =o X P \ Xr X PJ X P 


(A 17b) 
(A17c) 
(A 1 7d) 


My a = My + My X’ + My = 7 — + My — + (My + My \f) + My 9 (A17e) 

Y>A x t x u Y d v oo a v oo \ q « 


V 1 V T • / v . 

M17 A = M»7 , TT + lV^r 7 =-z — 4* tj (b + /IVI17 ” I\^r7 . “ ^/[.ry \b 

Z,A Z^ Voo Z^V m z p ^ z r z /3y Z /3 


(A17f) 


For small perturbation angles of attack a, the longitudinal forces can be expressed 
in terms of the aerodynamic lift and drag as 
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F x A = (La - D) (A 18a) 

F ZjA = -( L + Da ) (Ai8b) 

Taking the u, a, a, and q derivatives of equations (A18) and substituting these 
results into equations (A17a) and (A17c) gives the following expressions: 

f X,A - - D t - V - ( D o - L )\r - ( D a - L ) e < A19a > 

f Z,A - - L t - V - L d W- - ( L a ♦ D)# 1 - ( L q + L ci)8 - (La + ■>> (A 19b) 

v 00 v OO v ' 


where the higher-order perturbation terms are neglected. 

The aerodynamic forces and moments in equations (A19), (A17b), and (A17d) to 
(A17f) are nondimensionalized by the methods used in reference 9. The resulting aero- 
dynamic coefficients in these equations are then evaluated at trim conditions and the 
higher-order perturbation terms are neglected. The following relationships are thus 
obtained: 


X,A 


^( 2C D ♦ C D U ) 


X’ - 


pV S . _ , 

- °L 


z’ - 


pV Js 

^-( C La ' C L) 


pV 2 S 


f y,a - (^r c y^' + c y^V + pr 2 - c y p )^ 


PV to Sc 


— ( c ^ - cy P ) 


* 


Z,A 


^( 2C L * %) 


pSc 

~^L. 
4 a. 


pV S. 

-F^La * C d) 


Z’ 


(A20a) 




(A 20b) 


pV Sc, 


e - 


pV 2 s 

— 2~(° L a + Cd ) 


e - 


pvjs 


(A20c) 
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M 


X,A 


pS(c) 2 „ 
~~4 Cl r: 


pV Sc \ 

y + N^ c i> + 


PV ocS(g) Z c 


< P ~ 


pv^m 2 


(% - 


* 


lpV 2 Sc 


C,> 


(A20d) 


M- 


Y,A 


P V °° Sc he + c \ 

2 pm + L m u J 


x' + 


pS(c ) 2 „ 

4 '-m 


a 


z' + 


fpVc.Sc 


C m )z 


7 1 


at/ 


P V » S W 2 / C + c \ 

4 ( Cm a + C “q) 


0 + 


pV 2 Sc 
L—££ C 


m 


0! 


e + 


PV 0 q 2 Sc c 


m 


(A20e) 


M 


Z,A 


ps(e-) 2 „ 

> V °° SC C V + 

pVoo S(c) 2 „ 

_ 

pVooS(C) 2 ^ „ r 

4 n pJ 

y \ 2 v 

4 n P 


_ 4 (S ■ C ”r)J 






(A20f) 


Tether-Cable Forces and Moments 

The basic expressions for the tether-cable forces, given in reference 2, are derived 
in appendix B. In terms of the earth-fixed coordinate system, they are 


F X’,C = T 1 cos r 1 - k xx x c - k xz z c 

(A21a) 

f y',c = ' k yy y c 

(A21b) 

F Z',C = T 1 Sin y l " ^x x c " k zz z c 

(A21c) 

Assuming small perturbation angles, these forces can be written 
tem (fig. 3) as follows 

in the stability axis sys- 

F X,C = F X',C + f y',c^ " F z’,c 0 

(A22a) 

f y,c = " F X’,C^ + F Y',C + F Z',C 0 

(A22b) 

F z,c = F x',c 0 " f y’,c^ + F z',c 

(A22c) 
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In order to obtain moments about the balloon center of mass in the stability axis 
system, the coordinates of the bridle confluence point with respect to the center of mass 
must be determined. By inspection (fig. 4) these are 

\ = (*tr - *cg) cos a t + ( l tr - h cg) sin “ t 

(A23a) 

\ = (ttr - h cg) cos «t - [hr - Z cg) sin “t 
/ ^ 

(A23b) 

Therefore, the tether-cable moments about the balloon center of mass in the stability 
axis system are 

M X,C = ” h k 2 F Y,C 

(A24a) 

M Y,C = h k 2 F X,C " h k 1 F Z,C 

(A24b) 

M Z,C = h k 1 F Y,C 

(A24c) 

The displacements xj,, y^, and z' c of the upper end of the cable can be expressed 
in terms of the displacement of the balloon center of mass plus the rotation of the bridle 
confluence point about the center of mass. For small perturbations these displacements 
are 

x c = x ' + h k/ 

(A 2 5a) 

y c = y ’ + 

(A25b) 

z c = z ' - h k x e 

(A25c) 


Substituting equations (A21) and (A25) into equations (A22) and (A24) and neglecting 
higher order terms gives 

F X, C = ' k xx*' “ k xz z ’ - ( k x0 + T i sin + T ! cos Y x (A 2 6a) 

F Y, C = " k yy y ' + ( T 1 sin r i ‘ ^0)^ ■ ( T 1 cos y \ + ^0)^ (A 2 6b) 

F Z, C = * k zx x ' - k zz z ’ + ( T i cos 7^ - k zQ y + Tj sin (A26c) 
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where 


M x,c = - V - ( h k 2 T i sin y x + + ( h k 2 T i cos 

y i • VK 

(A26d) 

m y,c = - k 0x x ’ - k 0z z ' - k ee e - S' 1 ! sin V x + h k 2 T i 

cos Vl 

(A26e) 

M Z ,C = - V y ’ + ( h kj T l sin Yi - %$)* - (bk^i cos 

n + k w)' t ' 

(A26f) 

k x0 = h k 2 k xx ' h kj k xz 


(A27a) 

k z6> = h k 2 k zx " h kj k zz 


(A27b) 

k 0x = h k 2 k xx " h kj k zx 


(A27c) 

k 0z = fys^^xz ' k kj k zz 


(A27d) 

k ee = k 00 D + k ee T 


(A27e) 

k 00 D = h k 2 kxx ” h k 2 h k 1 ( k xz + k zx) + hj^zz 


(A27f) 

k 00 T = h k 2 ( T l Sin y l) + h ki( T l cos y i) 


(A27g) 

*y<l> = “bkg^y 


(A27h) 

k y 4> = bkj^y 


(A27i) 

k <*>y = 


(A27j) 

k o<f ’ h ' 2 k yy 


(A27k) 

k (f)x }/ = “ h k 1 h k 2 k yy 


(A271) 

v = 


(A27m) 

k 4/<P ~ k (p\ p 


(A27n) 

k M = bk^y 


(A27o) 


34 



APPENDIX A - Continued 


Buoyancy Forces and Moments 

The expressions for the buoyancy forces and moments about the center of mass in 
the stability axis system can be determined by inspection (figs. 1,3, and 4). Assuming 
small perturbation angles, these are 


F X,B " Be 


(A28a) 

F Y,B " ~ Bc f ) 

(A 2 8b) 

f z,b = -B 


(A28c) 

m x,b = - B 

( h cg - h br) cos “t + (^br “ Z cg) sin “t]* 

(A28d) 

m y,b = B ( 

\r ~ l cg) cos <*t - ( h cg - h br) sm “t] 


- B [( h cg - h br) cos a t + (^br - Z cg) sin “t] 0 

(A28e) 

m z,b = “ B 

(*br - Z cg) cos a t - ( h cg ' h br) sin “tj* 

(A28f) 


Gravity Forces and Moments 

In deriving the gravity forces for the equations of motion, only the components due 
to the structural weight of the balloon need to be considered. The effects of the apparent 
mass and the lifting gas are already included in the coefficients of the acceleration and 
buoyancy terms, respectively. 

The forces and moments due to gravity are determined by inspection (figs. 1, 3, and 
4) and for small perturbation angles are 


F X,G " " W s e 

(A29a) 

F Y,G = W s^ 

(A 2 9b) 

F Z,G = W s 

(A 2 9c) 

M X,G = - w s|( h sr - h cg) cos a t + (^sr + *cg) sin “t]^ 

(A29d) 


35 



APPENDIX A - Continued 


m Y,G = W s[0 

-sr + *cg) cos “ t ‘ ( h sr - h cg) sin “t] 


- w s 

[( h sr - h cg) cos “t + (*sr + z cg) sin “t] 0 

(A29e) 

m Z,G= - W s[ 

(^sr + ^cg) cos ®t ” (^sr “ ^cg) s ^ n 

(A29f) 


Balloon Equations of Motion 



If the expressions for each of the external forces and moments (eqs. (A20), (A26), 
(A28), and (A29)) are substituted into equations (A14) and (A15), the resulting equations of 
motion about the balloon center of mass are as follows: 

x-force 



+ B - W s - Tj sin y = 0 (A30c) 
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Rolling moment 


ps(tr) 


4 


y' - 


PV^Sc 


~2 If 


r + ( k ^ y )y + ( J x)^ - 

pv^so 2 


PVooS (tr)2^ ' 

3 % 


4> 


\ T l sin y x + k <^ + M s^> - M* + 


A | C 7 - C; 

4 ( ^ ‘r 


* 


IpVjSc 


2 C lp~\ T l C ° S r l + k ^ = 


= 0 


Pitching moment 


— ( 2C m + c m u ) x ' + ( k 0 X ) x ' " 


pS(c^ 

4 L m 


a 


■-rf - c »> + w 


W 


PV«.S(e) 2 fc ' 

'-m . + '-mqj 


' m d 


PV do 2 Sc 


6 + (k 00 + M Si 1 c m Q! 1^ 


PVoo 2 Sc c 


+ h ki Tj sin ?1 - h^Tj cos ?1 - M S2 = 0 


Yawing moment 


m 


pS(c ) 2 „ 

4 Cn /3 


y' ■ r“~f — c nJy' + Otyy) y ' ' ( Ixz )^ ■ 


pVooS(c) 2 


pv M s(c-) 2 „ 

4 ° n f 


0 


+ (m S2 + - hj^Tj sin y^cp + (I z )i£ + ^ (C n ^ - C nr j 


* 


fpvJSc 


+ I C n + h ki T 1 cos y x + = 0 


where 


- (^r ” z cg) B + (^sr + ^cg) w sj sin a t + |( k cg ” k br) B + ( h sr ' h cg) w sJ C0S a t 


M s 2 = (^br ' *cg) B + ( l sr + *cg) w s cos « t “ [( h cg " h br) B + ( h sr " h cg) w s] sin a t 


(A30d) 

(A30e) 


(A 3 Of) 

(A3 la) 
(A3 lb) 
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and 


m x = m x,o 


(A32a) 


m y “ m y,o 



(A32b) 


m z 


= m 


z,o 


+ 


PSc 

T" C L. 

a 


(A32c) 
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TETHER-CABLE CONFIGURATION AND FORCES 

The analysis of the cable shape and the forces exerted by the cable on the balloon 
bridle given in appendix B is adapted from the analysis of S. Neumark (ref. 2). The 
coordinate system and the forces acting on the cable are shown m the following sketch: 



Note that this coordinate system of appendix B is different from those employed elsewhere 
in the paper. 

Close examination of the expressions given in appendix B will show that many of the 
equations are invalid if either the cable weight per unit length w c or cable drag per umt 
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length n is zero. The proper expressions in these limiting cases are given in 
reference 2. 


Equilibrium Cable Shape 

The forces acting on an element of the cable of length dZ (see sketch) are the ten- 
sion T, cable weight w c d l, and cable drag P n dl normal to the cable. Drag along the 
cable is neglected in this analysis. The normal drag force per unit length P n depends 
on the component of wind velocity normal to the cable V n , the cable drag coefficient 
C D , and cable diameter d c as follows: 

Pn'CD^c lpV„2 (Bl) 

Retaining only first-order terms in the infinitesimals, the drag on the element is 

P n dl = n d l sin 2 y (B2) 

where 

n = C D c d c|pV 00 2 (B3) 

is the cable drag per unit length for a cable normal to the wind. Summing the x- and 
z-forces and retaining only first-order terms yields 

' , t 

P n dj sin y - T dy sin y + dT cos y = 0 (B4a) 

-w c dZ - P n d l cos y + T dy cos y + dT sin y = 0 (B4b) 

These equations may be combined to give 

dT = w c dZ sin y (B5) 

T dy = (n sin 2 y + w c cos y^dZ (B6) 

Dividing these equations gives 


dT 

T 


Sin 2 y + 2p cc 


(B7) 


where p = — Using the substitution f = cos y and partial-fraction decomposition in 
2n 

equation (B7) to obtain 
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dT _ p / df df \ 

T "q^q + p-f q-p + fj 


where q = yl + (p)^, and integrating from the upper end of the cable yields 

_1 

T 


(B8) 


(B9a) 


where 


r(y) = 


q + p - cos y 
q - p + cos y. 


P/q 


(B9b) 


For present purposes it is assumed that the cable length l, drag per unit length n, 

weight per unit length w c , tension at the upper end Tj, and angle of the upper end y 

are known. The procedure required to determine the coordinates of the upper end of 

the cable Xj and z^, the tension at the lower end Tq, and the angle of the lower end 

y is outlined as follows : Substituting equations (B9) m equation (B6) gives 
0 • 

T 

cU = -i i_ d y (BIO) 

nT l sin^y + 2p cos y | 

1 

Integration of this equation from the lower end to upper end yields 


l 




(Bll) 


where 


A(y) = r - d y . (B12) 

J 0 sin^y + 2p cos y 

X Q = X^yQj> and The angle y^ is unknown, but may be obtained by solving 

equation (Bll) for 

X o = x i - (B13) 

and using this value with equation (B12) to obtain 
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-j; 


T(y) 


® Jq sin^y + 2p COS y 


dy 


(B14) 


This equation is solved for the unknown limit of integration y^ by Newton iteration. 
With y^ known, equation (B9) is used to find 


rp _ rp 0 
L 0 ~ A 1 


(B15) 


where t q = r (^o)' From the sketch of appendix B, dx = dZ cos y. Using equation (BIO) 
in this expression yields 


dx = da 

nT l 


(B16) 


where da = 


r cos y 


sin^y + 2p cos y 


dy. Equation (B16) may be integrated numerically to give 


ri 

1 nT l J v„ 


da 


(B17) 


Finally, from the sketch and equation (B5) 


~ HT 

dz = sin y = 

w c 


(B18) 


which is integrated to give 


Zi = 


T i- T o 

w„ 


(B19) 


Cable Force Derivatives 

Consider the cable in its equilibrium position. If the upper end is slowly displaced 
in the XZ-plane from its original position (xj,z^ to a new position (xj + dx, Zj + dz), the 
resultant x- and z -force increments are 

dF x = k xx dx + k xz dz (B20a) 

dF Z - k zx dx + k zz dz (B20b) 
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where kxx, k xz , k zx , and k zz are the cable derivatives (spring constants) for the 
longitudinal case. The derivation of the expressions for these derivatives is given in 
reference 2. The results for the derivatives are 


‘'■xz 


= ±|Ti COS y^c 


k = — 
K zz - 6 


- sin y Q ) + n(z x - l sin y^sinSy^J 


(B21a) 

- cos y^J + n (l cos y^ - x^sin3y^ 

- 

(B21b) 

- sin - (w c + n sinSy^ cos y^J( 

- 1 sin ^ 

(B21c) 

- cos y^ - ^w c + n sin2y^ cos y^ 

{ l cos y 0 - *0 

(B21d) 


where 


6 = x^sin 7l - sin y Q j + z^cos y Q - cos yj - l sin^ - y Q j 


(B22) 


These derivatives depend only on quantities which were determined in the first section of 
this appendix. 

The single lateral cable derivative is determined by considering a small force dFy 
to act in the y-direction on the upper end of the cable. This will result in a small 
y-deflection according to 


dFy = kyydy (B23) 

where kyy is the cable derivative (spring constant) for the lateral case. The expression 
for this derivative is derived in reference 2 and is 


k yy " 


The value of kyy 


Ti( 

sin^y^ + 2p cos y^j 

C\ 

T (y) 

sin^y + 2p cos y 


is determined from equation (B24) by numerical integration. 


(B24) 
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TRANSFER OF AERODYNAMIC COEFFICIENTS WITH 
DISPLACEMENT OF AXIS- SYSTEM ORIGIN 

The equations of motion (eqs. (A30)) are written about the balloon center of mass 
which requires that the aerodynamic coefficients be referenced to this point. As the cen- 
ter of mass may change for various conditions, the aerodynamic coefficients are more 
conveniently specified about a fixed reference point and transferred to the mass center 
for each specific case. The following relations are developed for transferring the aero- 
dynamic coefficients from the arbitrary reference point to the center of mass. From 
figure 4 the geometrical relationships are 

x t = * C g cos a t + h cg sin a t 
z t = h cg cos a t - Z cg sin a t 

where x t and z t are the distances from the reference point to the center of mass. 

This development is given for completeness and to delineate some of the assump- 
tions involved. The quantities pertaining to the reference point are indicated by the sub- 
script R, and quantities pertaining to the center of mass are without subscripts. Only 
stability derivatives involving rates and displacements are considered. The apparent 
masses and inertias which involve accelerations are considered with the structural mass 
and inertias (e.g., eqs. (A10)). 


Longitudinal Coefficients 

The force and moment coefficients about the center of mass are determined by con- 
sidering the motion at the reference point resulting from motion about the center of mass 
and then calculating the resulting forces and moments at the center of mass due to the 
forces and moments acting at the reference point. 

Motion transfer .- The perturbations at the reference point generated by perturba- 
tions about the center of mass are 


R= 0 

(Cla) 

R = ® 

(Clb) 

R = ® 

(Clc) 


44 



APPENDIX C - Continued 


Z^ = Z' + X t 0 

(Cld) 

Z^ = Z’ + X t $ 

(Cle) 

= X' - Z t 6 

(Clf) 

Z' + Xi.$ 

a R =e+ V 

* oo 

(Clg) 

. Z f + Xfd 

“R= e+ V 

(Clh) 


Force and moment transfer.- The forces and moments at the center of mass in 


terms of those at the reference point are 

F Z = F Z,R ( C2a ^ 

F X = F X,R ( C2b ^ 

M y = m y,R ' z t M Z,R + x t M Z,R ^ C2c ^ 

Coefficients of lift.- For use in the z-equation the coefficients for 

( f z,a) r = -( L + d “r) < c3 > 


are needed. The procedure followed is: first, write equation (A17c) about the reference 
point; then, substitute equations (Cl) into the results; and finally, use the results of the 
substitution in equations (C2a) and (C3). Neglecting higher order terms, the resulting 
equation is equivalent to equation (A17c) but contains only independent variables per- 
taining to the center of mass and coefficients pertaining to the reference point. By 
equating coefficients of the independent variables of the resulting expression to the 
corresponding ones of equation (A 17c) and by evaluating the results at the trim condi- 
tion, the following transfer relations result: 


L 0! CL C¥,R 

(C4a) 

L . = C L . 

(C4b) 

a a,R 


L u " ° L u,R 

(C4c) 
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APPENDIX C - Continued 


C L„ " C L 


+ zxt c _fft/2C L + C L \ 
q,R c L a c l L ^u) 


Thus, the only term affected by the transfer is C T . 

q 

Coefficients of drag.- For use in the x-equation the coefficients for 


( f x,a), 


= Lo! t, - D 


are needed. Using equations (Cl), (C2b), (C5), and (A17a) gives 

C D = C D,R ( ( 

Cd « = Cd o,R (< 

Cd u = C D u ,r (< 

Thus, the aerodynamic drag coefficients are unaffected by equations (C6). 

Coefficients of pitching moment .- Using equations (Cl), (C2c), and (A17e) gives 


C m =C m,R + ^ C D ' = C L 

z t x t 

Cm u = C ”u,R*F C Ou'F Cl u 

' C m a> R - T( C L - C D„) - L„ * C D) 


c m • = c m • R _ ■=■ C L . 
a q!,R c a 


Cm q " Cm q,R " ^( CL q,R " 2Cm a) ' "r( 2Cm + Cm u) 


where terms such as C Y are neglected. 

^q 

Lateral Coefficients 

Motion transfer .- For lateral perturbations about the center of mass, the corre- 
sponding perturbations at the reference point are 
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APPENDIX C - Continued 


y^ = y' + : 

Z t 0 

- x t^ 


(C8a) 

H 

Ph 

•©- 




(C8b) 

4> r = 4> 




(C8c) 

II 




(C8d) 

"S>- 

II 




(C8e) 

(3 r =-^ + 

y’ 

Voo 

z t 4> 

+ V 

oo 

x t ^ 

’ V 

OO 

(C8f) 

(3 r = -i// + 

y' 

z t 4> 

+ V„ ' 

xt « P 

(C8g) 


Force and moment transfer .- The forces and moments at the center of mass in 
terms of those at the reference point are 


F Y = f Y,R 

M Z = M Z,R ‘ x t F Y,R 

M X = M X,R + z t F Y,R 

Coefficients of y-force .- Using equations (C8), (C9a), and (A17b) gives 


(C9a) 

(C9b) 

(C9c) 



(ClOa) 
(Cl Ob) 

(ClOc) 

(ClOd) 


Thus, only C v and C v are affected by the transfer equations (CIO). 

*p x r 

Coefficients of yawing moment.- Using equations (C8), (C9b), and (A17f) gives 
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APPENDIX C - Concluded 



l 0,R ' 1 % 

(Cl la) 

L - — C Y 

$,R c Y 0 

(Cl lb) 

'r.R' ffY r R +2 %) 

(Cllc) 

Xt 2z* 

- — Cv + — C n 

p,R c Y p,R c n /3 

(Cl Id) 


In the development of equations (Cll) the transfers involving C n . are neglected except 

for C n± itself. ^ 

0 

Coefficients of rolling moment.- Using equations (C8), (C9c), and (A17d) gives 


z t 

Cl - Cl + — C y 
%R c *0 


Ci “ Cl + — C y 
p,R c Y /3 


C, = c, + §/C Y + 2C, \ 
t?,r c \ y p,R l pj 


2x t 

C, =C ; +^C Y --s^C, 
fc r L r,R c 4 r,R c 


In the development of equations (C12) the transfers involving Cj 


for C 7 itself. 

L P 


P 


(C12a) 

(C12b) 

(C12c) 

(C12d) 


are neglected except 
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TABLE I - PROPERTIES OF 7 64 -METER BALLOON IN REFERENCE CONFIGURATION 


[See fig 4] 


Aerostatic properties 

Distance from balloon nose to reference point, a, m 
Atmospheric density, p, kg/m3 

Component of distance from reference point to bridle confluence point, ^ r , m 

Component of distance from reference point to bridle confluence point, tt r , m 

Buoyancy force, B, N 

Component of distance from reference point to center of buoyancy, £jj r , m 

Component of distance from reference point to center of buoyancy, hjj r , m 

Reference area, S, m2 
Reference length, c, m 

Component of distance from reference point to center of mass of balloon structure, Z gr , m 

Component of distance from reference point to center of mass of balloon structure, h sr , m 

Structural weight, W s , N 
Volume, Vjj, m3 

Dynamic masses and moments of inertia 
Apparent air mass at = 0), m xa , kg 
Apparent air mass at ( a ^ = 0), m>y a , kg 
Apparent air mass at = Oj, m z a , kg 

Component of distance from reference point to center of mass (includes apparent mass), l c g, m 

Component of distance from reference point to center of mass (includes apparent mass), h C g, m 

Mass of balloon structure and inflation gas m T kg 
Rolling moment of inertia (at a t = 0), Ijoc, kg-m2 
Pitching moment of inertia, Iyy. kg-m2 
Yawing moment of inertia (at = 0Y I 2Z , kg-m2 
Tether-cable parameters 
Diameter, d c , m 
Length, l, m 

Weight per unit length, w c , N/m 
Normal drag coefficient, C Dc 

Aerodynamic derivatives (specified about force reference point) 


■ 0 82(a t - 0 023) - 5 02(a t - 0 023) 3 + 111 4(a, - 0 023)5 


U L 


C, = 0 82 - 15 06 (a^ - 0 023) 2 + 557 0(a t ■ 
0 089 
c LqR = 0 685 

CL U = ° 

C D = 0 0487 + 186 2(a t - 0 023)6 
C D =■ 1117 2(a t - 0 023)5 
Cn =0 


0 023)4 


'•'D 

*^u 

C m,R 


= -0 0106 + 0 1435 a 


c m R = 0 H35 

QT,rl 

C m = -0 026 
m a,R 

C m = -0 189 

C m =0 
m u,R 

C Yn = -0 82 

C v = -0 089 

C v „ = 0 4 94 sin a* 
Y p,R 1 


C Y =0 685 
Y r,R 


' -Y r 

Cj = -0 1435 sm 

P,R 

C, =0 

%n 

C. = -0 0237 

l p,K 

C» = -0 178 sm a* 
*r,R 1 

S,R = ‘° 1435 

S,R = 0026 

c n = -0 0641 sm 2ar t 
P.R 1 

C n R =-0 189 


4 72 
1 225 

3 44 
3 82 
190 
2 15 
0 

7 04 
7 64 
-0 66 

0 38 
108 
19 0 

5 11 
23 9 
23 9 

1 10 
0 11 
14 2 
16 1 
171 
164 

0 0141 
61 0 
0 343 
1 17 
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TABLE D3.- PARAMETERS CONSIDERED IN LONGITUDINAL-STABILITY 

TREND STUDY (FIGS. 14 TO 27) 


mrp 

*yy 

l cg 

h cg 

m x,a 

m z,a 

B 

*br 

I 

I, Xj, and Zj 

he 

hr and a t * 

d c 

d c and Xj . 

d c and z^ . 

w c 

^xx 

^XZ 

kzx 

^ZZ 

AC D 

ac l 

ACm 

A % • • ■ • 

A ° m a ' ■ ' 

„ • • • • 

a,R 

C m .... 

m q,R 


Figure 

14 

15 
16(a) 
16(b) 
17(a) 
17(b) 
18(a) 
18(b) 
19(a) 
19(b) 
20(a) 
20(b) 
21(a) 
21(b) 
21(c) 

22 

23(a) 

23(b) 

23(c) 

23(d) 

24(a) 

24(b) 

24(c) 

25(a) 

25(b) 

26 

27 
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TABLE IV.- LIST OF PARAMETERS CONSIDERED IN LATERAL -STABILITY 

TREND STUDY (FIGS. 36 TO 44) 


*zz • 
l cg • 
m y,a 
B . 

W • 

l 

ac d 

ac l 

AC V 

Y |3 

AC n 

n /3 


Figure 

36 

37 

38 
39(a) 
39(b) 

40 

41(a) 

41(b) 

42 

43 

44 
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Balloon center of mass 
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Figure 1.- Tethered-balloon coordinate systems and forces considered in stability analysis. 






Center of bouyancy 
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- Geometry of balloon 












, rad/sec 



(a) u) and 77 versus V^. 

Figure 6.- Variation of longitudinal stability characteristics with velocity for 
reference configuration (tests 12 and 15). 
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n. 1/sec 


(b) Root locus plot with velocity as a parameter. 
Figure 6.- Concluded. 
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is shown for each velocity indicated.) 





Figure 7.- Concluded. 


UU-LLLUjJjjLU Ill 1 1 1 III I I 1 1 I I I I 1 1 1 1 1 1 1 I I I I 1 1 




V , m/sec 

CD 

Figure 8.- Variation of longitudinal stability characteristics u> 
and 7 ] with for configuration of test 10. 
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. rad/ sec 



V^. m/sec 


Figure 9.- Variation of longitudinal stability characteristics c o 
and 7j 'with for configuration of test 11. 
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V^. m/sec 



0 4 8 12 16 20 24 28 32 36 40 44 48 52 


V^. m/sec 

Figure 10.- Variation of longitudinal stability characteristics u> 
and 7j with V w for configuration of test 13. 
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n. l/sec a, rad/sec 



0 4 8 12 1 6 20 24 28 32 36 40 44 48 52 


V m , m/sec 



0 4 8 12 16 20 24 28 32 36 40 44 48 52 


m/sec 

Figure 11.- Variation of longitudinal stability characteristics w 
and tj with for configuration of tests 14 and 18. 
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r|. 1/sec o ) , rad/sec 



0 4 8 12 16 20 24 28 32 36 40 44 48 52 



0 4 8 12 16 20 24 28 32 36 4 0 44 48 52 


V^, mfsec 

Figure 12.- Variation of longitudinal stability characteristics o> 
and 77 with V w for configuration of test 16. 
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n ' 1/sec a. rad /sec 



0 4 8 12 16 20 24 28 32 36 40 44 48 52 


V , m/sec 

CD 

Figure 13.- Variation of longitudinal stability characteristics cj 
and Tj with V x for configuration of test 17. 
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. m/ sec (o , rad/sec 
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Figure 16.- Effect of center-of-mass location on calculated longitudinal stability boundaries 




(a) Apparent mass m X;a . (b) Apparent mass m z>£ 

Figure 17.- Effect of apparent mass on calculated longitudinal stability boundaries 












Calculat i ons 
Test 15 (unstable) 



(a) Stability boundary . ( 5 ) Balloon position. 

Figure 19.- Effect of tether-cable length on calculated balloon position and on 
calculated longitudinal stability boundary. 









Calculat i ons 
Test 15 (unstable) 

Test 16 (marginally stable) 
Test 17 (marginally unstable) 


E 



aas/pej -cn 


77 


Figure 20.- Effect of forward location of on calculated longitudinal stability boundaries 

and on calculated trim angle of attack. 






Figure 21.- Concluded. 
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Figure 23.- Effect of change in each tether-cable derivative on calculated longitudinal stability boundaries. 




oas/pej ‘ co 


D3S/UI ' "A 
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33S/PBJ ‘(0 


D3S/UJ 


(c) Fraction of calculated derivative k zx . (d) Fraction of calculated derivative 

Figure 23.- Concluded. 
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(a) Change in drag coefficient. (b) Change in lift coefficient. 

Figure 24.- Effect of static aerodynamic coefficients on calculated longitudinal stability boundaries. 







Reference 




H. l/sec a , rad/ sec 



0 4 8 12 16 20 24 28 32 36 40 44 48 52 

V^, m/sec 

(a) a; and 77 versus V,^. 

Figure 28.- Variation of lateral stability characteristics with velocity for 
reference configuration (tests 12 and 15). 
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Figure 28.- Concluded. 
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V 



Mode la at V = 10 m/s; -.681. 


Mode lb at V = 10 m/s; X - -2.44. 



Mode la at V = 20 m/s; X- -1.23. 



Mode lb at V = 20 m/s; X - -4.95. 




Mode la at V = 30 m/s; A = -1.80. 


Mode lb at V = 30 m/s; A = -7.47. 
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(a) Modes la and lb. 

Figure 29.- Lateral modes of motion for reference configuration. 



Figure 29.- Continued 





Figure 29.- Concluded 






n. ^ sec a), rad/sec 





Figure 30.- Variation of lateral stability characteristics a> 
and 7] with for configuration of test 10. 
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r|. 1/sec u , rad/sec 




V^, m/sec 

Figure 31.- Variation of lateral stability characteristics u> 
and 77 with V ^ for configuration of test 11. 
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n. 1/ sec w. rad/sec 



0 4 8 12 16 20 24 28 32 36 40 44 48 52 

V^. m/sec 

Figure 32.- Variation of lateral stability characteristics o> 
and ij with for configuration of test 13. 
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n. l^sec w, rad/sec 
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Figure 33.- Variation of lateral stability characteristics oj and 
77 with for configuration of test 14. 



33S/PBJ ‘rn D8S/1 'Ll 
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Figure 34.- Variation of lateral stability characteristics a> 
and tj with V,*, for configuration of test 16. 
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n. Ikec w, rad/sec 



Figure 35.- Variation of lateral stability characteristics u> 
and 7 ] with V w for configuration of test 17. 
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, m/sec u), rad/sec 
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b.n br.m 

(a) Buoyancy force. (b) Location of buoyancy force. 

Figure 39.- Effect of buoyancy force on calculated lateral stability boundary. 



m/sec to t rad/ sec 


1.2 


.8 


.4 


0 I ! I L 




Figure 40.- Effect of tether- cable length on calculated lateral stability boundary. 
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(a) Change in C D . (b) Change in C L . 

Figure 41.- Effect of static aerodynamic coefficients on calculated lateral stability boundary 
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Figure 42.- Effect of 



on calculated lateral stability boundary 



= -0.82 for reference configuration 
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Figure 43.- Effects of including various degrees of freedom and of AC n 


on calculated lateral stability boundary j 'C n ^ = -0.0254 for reference 
configuration^. 


(3 
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A motion- picture film supplement L-1118 is available on loan. Requests will be 
filled in the order received. You will be notified of the approximate date scheduled. 

The film (16 mm, 9 min, black and white, silent) shows computer generated longi- 
tudinal and lateral modes of motion for the 7.64-meter balloon tethered in steady winds 
of various velocities. 

Film supplement L-1118 is available on request to: 

NASA Langley Research Center 

Att: Photographic Branch, Mail Stop 171 

Hampton, Va. 23665 


CUT_ 

Date 

Please send, on loan, copy of film supplement L-1118 to 

TN D-7272. 


Name of organization 
Street number 


City and State Zip code 

Attention: Mr. 

Title 


CUT 


Place 

Stamp 

Here 


NASA Langley Research Center 

Att: Photographic Branch, Mail Stop 171 

Hampton, Va. 23665 
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